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Abstract 

We consider the computation of the volume of the union of high-dimensional 
geometric objects. While showing that this problem is #P-hard already for 
very simple bodies, we give a fast FPRAS for all objects where one can (1) test 
whether a given point lies inside the object, (2) sample a point uniformly, and 
(3) calculate the volume of the object in polynomial time. It suffices to be able 
to answer all three questions approximately. We show that this holds for a large 
class of objects. It implies that Klee's measure problem can be approximated 
efficiently even though it is #P-hard and hence cannot be solved exactly in time 
polynomial in the number of dimensions unless P — NP. Our algorithm also 
allows to efHciently approximate the volume of the union of convex bodies given 
by weak membership oracles. 

For the analogous problem of the intersection of high-dimensional geometric 
objects we prove #P-hardness for boxes and show that there is no multiplicative 
polynomial-time 2^* -approximation for certain boxes unless NP — BPP, but 
give a simple additive polynomial-time e-approximation. 

1 Introduction 

Given n bodies in the c?-dimcnsional space, how efSciently can we compute the volume 
of the union and the intersection? We consider this basic geometric problem for 
different kinds of bodies. The tractability of this problem highly depends on the 
representation and the complexity of the given objects. For many classes of objects 
already computing the volume of one body can be hard. For example, calculating 
the volume of a polytope given either as a list of vertices or as a list of facets is 
#P-hard pin. For convex bodies given by a membership oracle one can also show 
that even though there can be no deterministic (0(1)^/ log c?)'*-approximation for 
> 2 Q , one can still approximate the volume by an FPRAS (fully polynomial-time 
randomized approximation scheme). In a seminal paper Dyer, Frieze, and Kannan [12 1 
gave an 0*{cP^) algorithm, which was subsequently improved in a series of papers [2, 
3, 2l|,[22| to [i^l (where the asterisk hides powers of the approximation ratio 



and logd). 

Volume computation of unions can be hard not only for bodies whose volume is 
hard to calculate. One famous example for this is Klee's Measure Problem (KMP). 
Given n axis-parallel boxes in the c?-dimensional space (d constant), the problem asks 
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for the measure of their union. In 1977, Victor Klee showed that it can be solved 
in time O(nlogn) for d = 1 This was generahzed to d > 1 dimensions by 

Bentley in the same year. He presented an algorithm which runs in 0{n'^~^ logn), 
which was later improved by van Leeuwen and Wood to 0{n'^~ ). In 1988, 
Overmars and Yap [IHj obtained an 0(n''/^logn) algorithm. This was the fastest 
al gorithm for d ^ 3 un til very recently Chan Q presented a slightly improved version 
of lOvermars and Yapl 's algorithm that runs in time n'^/^ 2'-'('°s where log* denotes 
the iterated logarithm. So far, the only known lower bound is J7(nlogn) for any d 13|. 
Chan also proves that no algorithm of runtime possible assuming W[l] ^ 

FPT, which is a weaker result than P NP and a commonly accepted conjecture on 
fixed-parameter tractability. Note that the worst-case combinatorial complexity (i.e., 
the number of faces of all dimensions on the boundary of the union) of Q{n'^) does 
not imply any bounds on the computational complexity. There are various algorithms 
for special cases, e.g., for hypercubes U, [l3| and unit hypercubes [1]. In this paper 
we explore the opposite direction and examine the union of more general geometric 
objects. 



Our results 

It is not hard to sec that KMP is #P-hard (sec Thcorcm[T]). Hence it cannot be solved 
in time polynomial in the number of dimensions unless P ~ NP. This shows that 
exact volume computation of unions is intractable for all classes of bodies that contain 
axis-parallel boxes. This motivates the development of approximation algorithms for 
the volume computation of unions. Based on an FPRAS for #DNF by Karp. Luby, 
and Madras [l8|, we give an FPRAS for a large class of bodies including boxes, 
spheres, polytopes, convex bodies determined by an oracle, and schlicht domains. 
Additionally, also fixed afSne transformations of the forementioned objects can be 
allowed. The underlying bodies B just have to support the following oracle queries 
in polynomial time: 

• PointQuery(x, B): Is point x € R'' an element of body B? 

• VOLUMEQuERY(i?): What is the volume of body B7 

• SampleQuery(S): Return a random uniformly distributed point x € B. 

PointQuery is a very natural condition which is fulfilled in almost all practical cases. 
The VolumeQuery condition is important as it could be the case that no efficient 
approximation of the volume of one of the bodies itself is possible. This, of course, 
prevents an efficient approximation of the union of such bodies. The SampleQuery 
is crucial for our FPRAS. In Section [2731 we will show that it is efficiently computable 
for a wide range of bodies. 

An important feature of our algorithm is that it suffices that all three oracles are 
weak. More precisely, we allow the following relaxation for every body B (vol(_B) 
denotes the volume of a body B in the standard Lebesgue measure on W^, more details 
are given in Section [5]): 

• PoiNTQuERY(a;, B) answers true if and only if x G B' for a fixed B' C M'^ with 
VOl{{B' \B)U{B\ B')) ^ epVOL(B). 

• VOLUMEQUERY(i?) returns a value V with (1 - ey) yol{B) ^ V ^ {I + 
ev) vol(B). 
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• SAMPLEQUERY(i3) returns only an almost uniformly distributed random 
point [16|, that is, it suffices to get a random point x € B' (with B' as 
above) such that for the probability density / we have for every point x: 
|/(a-)-l/vOL(B')| <es- 

Let P{d) be the worst PointQuery runtime^ of any of our bodies, analogously V{d) 
for VolumeQuery, and S{d) for SampleQuery. Then our FPRAS has a runtime 
of 0{nV{d) + {S{d) + P{d))) for producing an e-approximation^ with probability 

o 2 2 

> I if the errors of the underlying oracles arc small, i.e., es,sv < -^f^ and ep < -^f^- 
For example for boxes (e.g., for KMP), this reduces to 0(1?) and is the first FPRAS 
for this problems. In Section [231 we also show that our algorithm is an FPRAS for 
the computation of the volume of the union of convex bodies. 

The canonical next question is the computation of the volume of the intersection of 
bodies in R"*. It is clear that most of the problems from above apply to this question, 
too. #P-hardness for general, i.e., not necessarily axis-parallel, boxes follows directly 
from the hardness of computing the volume of a polytope [ill , [lij . This leaves open 
whether there are efficient approximation algorithms for the volume of intersection. In 
Section[3]we show that there cannot be a (deterministic or randomized) multiplicative 
2'^ -approximation in general, unless NP = BPP by identifying a hard subproblem. 
Instead we give an additive e-approximation, which is therefore the best we can hope 
for. It has a runtime of 0{nVid) + e'"^ Sid) + ne-'^P{d)), which gives for 
boxes. 



2 Volume computation of unions 

In this section we show that the volume computation of unions is #P-hard already 
for very simple axis-parallel boxes. After that we give an FPRAS for approximating 
the volume of the union of bodies which satisfy the three aforementioned oracles and 
describe several classes of objects for which the oracles can be answered efficiently. 



2.1 Computational complexity of union calculations 

Consider the following problem: Let 5 be a set of n axis-parallel boxes in R*^ of the 
form B = [ai, 6i] x • • • x [a^, bd] with a^, bi G R, < bi. The volume of one such box 
is V0L(i3) = niLi(^j ~ ^i)- To compute the volume of the union of these boxes is 
known as Klee's Measure Problem (KMP). 

It is know that the associated decision problem of deciding whether there is a point 
that is not in the union is NP-hard Q . We consider the actual counting problem and 
prove m the following Theorem[l]that KMP is #P-hard. To the best of our knowledge 
there is no published result that explicitly states that KMP is #P-hard. However, 



without mentioning this implication, Suzuki and Ibaraki j27| sketch a reduction from 
#SAT to KMP. We present a reduction from #MON-CNF to KMP which we can 
reuse in Theorem |6] for the hardness proof for intersections. #MON-CNF counts the 
number of satisfying assignments of a Boolean formula in conjunctive normal form in 



^The runtime of the oracles can also depend on the required approximation guarantees. In order 
to simplify the notation, this dependency is not made explicit. 
■^We will always assume that e is small, that is, < £ < 1. 
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which all variables are unnegated. While the problem of deciding satisfiability of such 
formula is trivial, counting the number of satisfying assignments is #P-liard pij . 

Theorem 1. KMP is #P-hard. 

Proof. To reduce #MON-CNF to KMP, let / = ALi VjeC;, be a monotone 
Boolean formula given in CNF with Ck C [d] := {1, . . . , d} , ioi k G [n], d the number 
of variables, n the number of clauses. Since the number of satisfying assignments 
of / is equal to 2'^ minus the number of satisfying assignments of its negation, we 
instead count the latter: Consider the negated formula / = Vfc=i Ai^Ck I^i'^^tj 
we construct a box Ak = [0,q[''^] x • • • x [0,q^J^^] in for each clause Cfc with one 
vertex at the origin and the opposite vertex at (qj'^'' , . . . , 9^'^''), where we set 

\ 2, if * ^ Cfe ' ^ ^ 

Observe that the union of the boxes can be written as a union of boxes of the form 
Ux = [xi,xi + 1] X • ■ • X [xci, Xd + 1] with X = (xi, . . . , x^) G {0, 1}''. Let x 6 {0, 1}'' 
and Ux C yJl^i Ak. Then there is a fc such that qp' = 2 for all i with Xi = 1. By 
definition of q^-'^'' , this implies that AiGCt ^^'^ ^^^^ f ^'"^ satisfied. 

The same holds in the opposite direction, that is, if x satisfies / then Ux Q 
Ufc=i Hence, since vol{Ux) = 1, we have vOL(U^^^Afc) = \{x e {0,1}'' | 
a;satisfics/} | . Thus a polynomial time algorithm for KMP would result in a polynomial 
time algorithm for ^^MON-CNF, which proves the claim. □ 

Note that we actually proved a little bit more than stated in the theorem. That 
is, we proved that even calculating the volume of the union of boxes which all have 
the point 0'' in common is ^^P-hard. This specific problem is known as hypervolume 
indicator [s^l and is a very popular and widely used measure of fitness of Pareto sets 
in evolutionary multi-objective optimization. 



2.2 Approximation algorithm for the volume of unions 

In this section we present an FPRAS for computing the volume of the union of 
objects for which we can answer PointQuery, VolumeQuery, and SampleQuery 
in polynomial time. The input of our algorithm ApproxUnion are the approximation 
ratio £ and the bodies Bi, . . . , Bn in M.'^ defined by the three oracles. It computes an 
approximation U S E of C/ := vol(1J"^-^ _Bi) such that 

Pr[(l--e)[/<J7<(l+e)[/] >| (1) 

in time polynomial in n, 1/e and the query runtimes. Note that the constant | can 
be increased to any number by using a probability amplification technique. 

We are following the algorithm of Karp et al. [18| which the authors used for 
approximating #DNF and other counting problems on discrete sets. The two main 
differences are that here we are handling continuous bodies in R'' and that we allow 
erroneous oracles. The latter relaxation is crucial to incorporate, amongst other 
things, the class of convex bodies. Our algorithm ApproxUnion is shown on page [3 
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Algorithm 1: ApproxUnion {S, e, ep, ey, es) calculates an e-approximation of 
U = vol(1J"^^ Bi) for a set of bodies S = {Bi, . . . , in M'' determined by 
the oracles PointQuery, VolumeQuery and SampleQuery with error ratios 
£p, £v, ^s- 



/V (l+£s) (l+gy) (1+nep) 

^ (l-£v)(l-Ep) 
rp 241n(2) (l+e)K 

■~ 1-2-8 (C-l)n 

for all i £ [n] do 

compute := VOLUMEQUERY(i3i) 

od 

V :=Er=iK' 

for M to oo do 

choose i € [n] with probability V- /V 

X := SAMPLEQUERY(_Bi) 

tM 
repeat 

if to + ... +tM >T then return ^ 

choose j £ [n] uniformly at random 

tM := + 1 
until PointQuery {x,Bj) 
od 



The total number of steps of the algorithm is T. This number is chosen in advance 
such that one can prove that equation ([1]) holds. The algorithm itself is very simple. 
First, it queries the volumes'^ of the bodies Bi and computes V = Ylll=i ^1 ■ Then 
it repeats the following: It chooses a body Bi with probability (roughly'') proportional 
to its volume and chooses a point x G Bi (roughly*) uniformly at random. Afterwards, 
the algorithm chooses bodies Bj with probability 1/n and (roughly*) checks whether 
X € Bj. The number tM of chosen bodies until we find a Bj with x G Bj can then be 
used to estimate how many bodies Bj cover x. 

For this, observe that the point x is chosen with probability density k{x)/V' with 
k{x) ~ \i £ [n] \ X £ Bi\. Each number tM has expected value (roughly*) n/k{x) 
for a fixed x G Ui^i- Hence when PointQuery {x,Bj) answers yes, tM should 
be of the order of (roughly*) n/k{x) ■ k{x)lV' dx = jJin/V) dx = nU/V. This 
implies that when U is returned and T = tg + . . . + tA/, the value nU M/V' is near 
T, i.e., TV' /nM is near U. This gives the intuition why the algorithm returns the 
approximation = 

In Section ^TM we show correctness of ApproxUnion. More precisely, we show 

that it returns an ^-approximation with probability > | and T = 0(72-) if £Sj£v 
22 

and ep < . The last inequality reflects the fact that we cannot be arbitrarily 
accurate if the given oracles are inaccurate. If all oracles can be calculated accurately, 

^Note that here and in the remainder an unprimed variable denotes an exact value and a primed 
variable denotes a value subject to some error introduced by the erroneous oracles, 
^"roughly" only for erroneous oracles. 
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i.e., if ep = es = £v = 0, the algorithm runs for just T = many steps. 

Overall, the runtime of ApproxUnion is clearly 0{n-V{d) + M-S{d) + T-P{d)) = 
0{n ■ V{d) + T ■ {S{d) + P{d))), where V{d) is the worst VolumeQuery time for 
any of the bodies, analogously S{d) for SampleQuery and P{d) for PointQuery. 
If es,ev < ^ and ep < the runtime is 0{nV{d) + ^ {S{d) + P{d))). 

For boxes all three oracles can be computed exactly in 0{d). This implies that 
our algorithm ApproxUnion gives an £-approximation of KMP with probability 
> I in runtime O(^). For more complex objects like convex bodies determined 
by a membership oracle, there are no exact oracles. The following section discusses 
different classes of objects for which our algorithm can be applied. 



2.3 Classes of objects supported by our FPRAS 

Finding an FPRAS for the union of a certain class of geometric objects now reduces 
to calculating the respective PointQuery, VolumeQuery and SampleQuery in 
polynomial time. We assume that we can get a random real number in constant 
time. Then all three oracles can be calculated in time 0{d) for d-dimensional boxes. 
This already yields an FPRAS for the volume of the union of arbitrary boxes, e.g., 
for KMP. Note that if we have a body for which we can answer all those queries, all 
affine transformations of this body fulfill these three oracles, too. We will now present 
three further classes of geometric objects. 



Generalized spheres and boxes Let be the class of boxes of dimension k, 
i.e., Bfc — {[ai,6i] x ■•• x [ak,hk] \ ai,bi 6 M, < bi} and the class of spheres 
of dimension i. We can combine any box B € and sphere 5* £ S^-fc to get a 
d-dimensional object B x S. Furthermore, we can permute the dimensions afterwards 
to get a generalized "box-sphere." in M.^ this corresponds to boxes, spheres and 
cylinders. VolumeQuery can be computed easily as we can compute the volume of 
a sphere by a well-known formula and thus the volume of the product B x S. As one 
can check whether a given point x = {xi, . . . ,Xd) lies in B x 5* by checking whether 
(xi, . . . ,Xk) lies in B and {xk+l^ ■ ■ ■ ,Xd) lies in S*, also PointQuery is a standard 
task of geometry. To answer SampleQuery, it suffices to choose a random point 
(xi, . . . ,Xk) in B and to choose a random point inside the sphere S', which can be 
done in polynomial time as described, e.g., by Muller [23 |. 



Convex bodies As mentioned in the introduction, exact calculation of Volume- 
Query for a polytope given as a list of vertices or facets is #P-hard [ll|, [3l ■ Since 
there are randomized approximation algorithms (see Dyer ct al. [12| for the first one) 
for the volume of a convex body determined by a membership oracle, we can an- 
swer VolumeQuery approximately. The same holds for SampleQuery as these 
algorithms make use of an almost uniform sampling method on convex bodies. See 
Lovasz and Vempala for a result showing that VolumeQuery can be answered 
with 0*(^) questions to the membership oracle and SampleQuery with O* (ir) 
queries, for arbitrary errors ev,£s > (where the asterisk hides factors of log(fi) 
and log(l/ev) or log(l/es))- PointQuery can naturally be answered with a single 

2 __ 

question to the membership oracle. By choosing ey = es = , Theorem [2] together 
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with Lemma [3] shows that ApproxUnion is an FPRAS for the volume of the union 
of convex bodies which uses membership queries. 

Star-shaped bodies Star-shaped bodies are a generahzation of convex bodies 
which have at least one point such that every line through the point has a convex in- 
tersection with the body. They can also be viewed as the union of convex sets, with all 
the convex sets having a nonempty intersection. The subset of points that can "see" 
the full set is called the kernel of the star-shaped set. Assuming that we are given 
membership oracles for the body as well as for the kernel, Chandrasckaran, Dadush, 
and Vempala [l^ recently showed that for star-shaped bodies SampleQuery can 
be answered with C'*(^^p-) questions to the membership oracle, where rj is the frac- 
tion of the volume taken up by the kernel. We can also approximate the volume of 
the (convex) kernel with 0*(ir) questions to the membership oracle as discussed 
above and estimate t} by O(-t-t) samples. Then the runtime of VolumeQuery 

is 0*(^ -\ — 5-t-t). PointQuery can again naturally be answered with a single 
question to the membership oracle. 

Schlicht domains Let continuous functions with at < hi , where 

oi, bi are constants. Let K dW^ he defined as the set of all points (xi, . . . , Xd) G 
such that ai < xi < &i,a2(a;i) < X2 < b2{xi), . . . , ad{xi, . . . , Xd-i) < Xd < 
bd{xi, . . . ,Xd-i)- K is called a schlicht domain in functional analysis. Fubini's the- 
orem for schlicht domains states that we can integrate a function /: — > R by 
iteratively integrating first over Xd, then over Xd-i-, until we reach xi. This 
way, by integrating /(■) — 1, we can compute the volume of a schlicht domain as 
long as the integrals are computable in polynomial time, and thus answer Volume- 
Query. Similarly, we can choose a random uniformly distributed point inside K: 
Let K{y) = {(xi, . . . ,Xd) & K \ xi — y}. Then K(y) is another schlicht domain for 
every oi < y < &i. Assume that we can determine the volume of every such K(jj) 
and the integral I{y) = J^^ K{x) dx. Then the inverse function : [0, V] — ?> R, 

where V = /^^^ K{x) dx is the volume of allows us to choose a y in [ai, hi] with 
probability proportional to vOL(_ft'(y)). By this we can iteratively choose a value y 
for xi and recurse to find a uniformly random point (7/2, ■ • ■ in K(y), plugging 
both together to get a uniformly distributed point (yi, . . . , yd) in K. Hence, as long as 
we can compute the involved integrals and inverse functions (or at least approximate 
them good enough), we can answer SampleQuery. Since PointQuery is trivially 
computable ~ as long as we can evaluate and hi efficiently - this gives an example 
showing that the classes of objects that fulfill our three conditions include not only 
convex bodies, but also certain schlicht domains. 

Note that all above mentioned classes of geometric objects contain boxes and hence 
our hardness results still hold and an e-approximation algorithm is the best one can 
hope for (unless P = NP). 



7 



2.4 Analysis of our algorithm 

We now show correctness of our algorithm ApproxUnion described in Section [ 
and prove bounds for its approximation ratio. The following theorem is our main 
result for ApproxUnion. It holds for exact and weak oracles. 

Theorem 2. Given errors < ep,es,ev < ^ of the queries, the algorithm 
ApproxUnion({_Bi, . . . , Bn}, £, £p, ev, es) returns a value U with 



choosing 



under the condition 



Pt[{1 - e)U <U < {l+e)U] > 

^_ 241n(2)(l + g)n 

£2 - 8 ((7 - 1) n ' 



e > 



ev + 2(l+ev)V^2(C'- (2) 



where U vOL(Ur=i B.) , e f^, and C := ^'^'^^^It^J^lltr' ' 

First note that with accurate oracles, i.e., if ep = £s = £v = 0, we get e = e, 
(7=1 and, thus, T = 24in(2) (i+e) n ^ condition ^ becomes trivial, above 

theorem implies that our algorithm is indeed an FPRAS. 

Given non-zero query errors, one clearly cannot be arbitrary accurate, which is 
reflected by the lower bound ([5]) for e. However, the following lemma shows that 
condition ([2|) is fulfilled for small enough ep, eg and ey- 

Lemma 3. Fores,£v < e^/(i7n) andep < e^/{4:7n'^) the condition ([2]) of Theorem\^ 
is fulfilled and we have T = ©(p-) . 

Proof. As the right hand side of ([2|) is increasing in ep, ey and eg, we can assume 
w.l.o.g. that ey^= = and ep = This gives C < (l + ^) (l - ^) . 

Observe that < jj. Since (1 + x)^ (1 - x)"^ < 1 + ax holds for a := and 

also < a; < 4^, we have C < 1 + Q.-^. Hence, we can upper bound the right hand 
side of @ by 



ey + 2(l + ey)v'2((7-l)n < ^ + 2 (l + ^) ^J2a^ 

^^+2(1+^) v¥) 



as 1/n < 1 and e < 1. Since ^ + 2(1 + < 1, condition ^ is fulfilled. 

We now bound the terms e and C. Since ey > 0, we clearly have e < e implying 

2 

1 + e < 2. Furthermore, since ey < it also holds that 



^ £- ^ ^ £(47n ~ e) ^ e46n _ 23 



l + £v~l + -i^ + 4771 - 48?! 24 

47n 
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where we used e < 1 and n > 1. Using this and the upper bound for C we get for the 
denominator of T: 



Therefore, we get 



c o(C — l)n 4575168^^ 

In order to prove Theorem [2l we will generalize the corresponding proof of Karp 
et al. to cover weak oracles. For most lemmas it suffices to insert the error 
constants ep, es, ev, but for the main proof of Theorem [2] one has to be little bit more 
careful. 

First, recall that we are given bodies Bi,...,Bn by oracles, where 
PoiNTQuERY(a;, returns true for every x £ B'^, such that the result is wrong 
for the set W., = {B^\B'^ U {B[\Bi) with MOi.{W{) < epVOh{B,), which implies 

(1 - ep) vol(B,) < vOL{Bl) < (1 + ep)vOh{B,). (3) 

Furthermore, the volume Vl of body Bi is computed by VolumeQuery. VJ is an 
ey-approximation of the corresponding exact volume Vi, i.e., 

il~ev)V,<V^ <{l+ey)V,. (4) 

We set V := ^"^i V/ and V J27=i Then it clearly holds that 

il-sv)V <V' <{l+eY)V. (5) 

Furthermore, let U be the exact volume of the union of the B^'s and /.i = U/V. 

As in Karp et al. [l^l, we define for a point x G R'' the number of covering bodies 
COV(a;) = |{i e [n] I PointQuery(2;, _B,) — truc}|. Additionally, we set 

Rk := {a; e R'' I COv(a;) = k} 

and rk := VOL{Rk). Then we have Ylk=i ^'''k = I]"=i vol(B-), so that 

n 

(l-ep)T^<^fcrfc < (l + £p)K (6) 

k=l 

Furthermore, J2k=i ''fc = ^'^^(\J7=i ^i)' ^^at 



(1 - nep) (7 < y rfe < (1 + nep) U. (7) 



fc=i 



To get a sample point in our algorithm, we first choose an i G [n] with probability 
VI /V and then choose a random point x in B[ via SAMPLEQuERY(_Bi). We consider 
the probability of this random point x to lie in the region Rk ■ With error- free oracles 
this probability would be Pr[a; G Rk] ~ ^r- as exactly k bodies cover each point of Rk 
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and we have ^''"k = V m the error-free setting. With errors, simple calculations 

using inequalities ([3]), (|4]) and ([5]) yield 

(l-£s)(l-£v) krk , jl + es) (1 +£v) krk 

(l + £v)(l + sp) — ^ P'^t" ^ ^ (l-£v)(l-ep) — ■ 

In the algorithm, denotes the number of iterations in the inner loop during the 
m-th iteration of the main loop, i.e., the number of trials to find a box containing 
the m-th point x. These tm are independent identically distributed random variables. 
Let i be a variable distributed as each tm- Then the following Lemma holds. 

Lemma 4. Lei < A < i and C ^ (i+y;) (i+^v) (i+"^p) . Then 

— —2' (l-ev)(l-ep) 

Ex[e^*/"] < (7e(^+^^')^ and 
Ex[e-^*/"] < Ce-^^-^">''. 

These bounds closely match Lemmas 7 and 9 of Karp et al. Adapting the 

proof is straightforward. The factor C arises by the usage of inequalities ([7]) and (|5]). 

In the remainder, let Se := X]i=o ^"-^ ^^'^ step at which the £-th trial is completed 
and let Ni be the number of trials completed after step i, so that in the end M = Nt- 
Then Ni < £ if and only if Se > i. 

Corollary 5. Lete <2. Then 

Vy[Si > (1 + £) n^it] < C'^e-""''^/^ and 
Pr[Se < (1 - e)n^i£] < C^e-^^'^/*^. 

This corresponds to Corollaries 8 and 10 in Karp et al. [l3|- We can reuse their 
proof word for word; the only change is our Lemma |4] which brings in the factor C. 

It remains to prove Theorem [2j The corresponding theorem of Karp et al. [3l is 
called Self- Adjusting Coverage Algorithm Theorem II. However, adapting their proof 
is not as straightforward as for the previous lemmas. It is presented in more detail in 
the remainder of this section. 

Proof of Theorem^ Let £ = and C = and assume 



£>£v + 2(H-£v)V2(C-l)n. (9) 

This implies £ > and £^ - 8 (C* - 1) n > 0. Let 

^ 241n(2)(l-F£) 
^ At(£2-8(C-l)n)(l + £) 

and 

241n(2) {l+e) 



k2 



Ai (£2 - 8 (C - 1) n) (1 - £) ' 
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Then T = fcin/i (1 + e) = fc2?^/i (1 — £)■ Now, if we have fci < M < k2, then 

T T T 
k2 ~ M ~ ki 

and thus 

TV TV ~ TV 
< < . 

nk2 nM nki 

By plugging in T, fci and fc2 we get 

il-s)^iV <U < {l+e)fiV. 

A httlc calculus shows that (1 — e) (1 — ev) > 1 — e, only based on the definition of 
£, the non-negativity of ey and e, and sy < s. By using this, equation ([S]), the fact 
(1 + e) (1 + ev) = 1 + £, and = U/V we get 

il-e)U <U <il + e)U, 

and thus the estimation is an e-approximation, if fci < M < k2. Hence, it suffices to 
show 

Pr[Af > fcs] + Pr[Af < fci] < \. (10) 

We have 

Pr[M < fci] = Pr[^fe, > T] = Pr[5fc, > k,n^l (1 + e)]. 
By £ > we have, using Corollary El 

Pr[Af < fci] < Cl'i^-f^S'>'i/S < 5(0-1) 

By inequality ([9]) we get /i£^ — 8(C' — l)n/i > and since /.t > 1/n we have iii^ — 
8 {C -1) > ^ie^ -8{C - 1) riAi- Hence, 

Ft[M < fci] < e-3i"(2) ^ 1^ 

We analogously get Pr[Af > fc2] < |- Plugging both results in equation ([TU]) finishes 
the proof. □ 



3 Volume computation of intersections 

In this section we are considering the complement to the union problem. We show 
that surprisingly the volume of a intersection of a set of bodies is often much harder 
to calculate than its union. For many classes of geometric objects there is even no 
randomized approximation possible. 

As the problem of computing the volume of a polytope is #P-hard [HI, [3] , so is 
the computation of the volume of the intersection of general (i.e., not necessarily axis- 
parallel) boxes in M.'^. This can be seen by describing a polytope as an intersection of 
halfplanes and representing these as general boxes. 
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Let us now consider the convex bodies again. Trivially, the intersection of convex 
bodies is convex itself. From the oracles defining the given bodies Bi,. . . ,Bn one can 
simply construct an oracle which answers PointQuery for the intersection of those 
objects: Given a point a; 6 M'' it asks all n oracles and returns true if and only if x lies 
in all the bodies. One could now believe that we can apply the result of Dyer et al. 



12| and the subsequent improvements mentioned in the introduction to approximate 
the volume of the intersection and get an FPRAS for the problem at hand. The 
problem with this approach is that the intersection is not "well-guaranteed." That is, 
there is no point known that lies in the intersection, not to speak of a sphere inside 
it. However, the algorithm of Dyer et al. [l2j relies vitally on the assumption that 
the given body is well-guaranteed and hence cannot be applied for approximating the 
volume of the intersection of convex bodies. 

We will now present a hard subproblem which shows that the volume of the 
intersection cannot be approximated (deterministic or randomized) in general. 

Definition 1. For p e [0, 1]'', let {x | ^ =^ ^ VI < i < d}. We call 

p 



Bp [0,1]'^ \ Bp a co-box 



A co-box is a box where we cut out another box at one corner. The resulting object 
can itself be a box, too, but in general it is not even convex. It can be seen as the 
complement of a box Bp relative to a larger background box [0,1]''. Note that it is easy 
to calculate the union of a set of co-boxes {Bp^, Bp^, . . . , i?p„} with pi,p2, ■ ■ ■ ,Pn G 
[0, 1]'' as U Bp. ~ [0, l]''\n -Bpi- On the other hand, the calculation of the intersection 
of a set of co-boxes is 

#P-hardbyThcorcm[T]asU5p. = [0, l]''\n ^p. ■ The following 
theorem shows that it is not even approximable. 

Theorem 6. Let pi,p2, ■ ■ ■ ,Pn £ ^>q- Then the volume of n"=i cannot be 
approximated in (deterministic or randomized) polynomial time by a factor of 2'^ 
for any e > unless NP = BPP. 

Proof. Consider again the problem #MON-CNF already defined in Section [21 Let 
/ = Afe=i y ii^Ck ■'^^ ^ monotone Boolean formula given in CNF as defined in 
the proof of Theorem [1] Wc now construct for every clause Ck a co-box Bp,, with 

Pk = {p'i\ ■ ■ ■ jpli^) E^nd p[''^ = 1/2 if i e Ck, and p^'^^ = 1 otherwise. The boxes Bp,^ 
here correspond to the Ak in the proof of Theorem [TJ For x = {xi,. . . , Xd) G {0, l}"^ 
let Ux = \{xi + \)\ X • • • X \_\xd, \{xd + 1)] • As in the proof of Theorem [1] one 

shows that C {Sl^^ Bp^ if and only if x satisfies / and hence also Ux C Hfe^i if 
and only if x satisfies /. This implies that the volume of nr=i ^Pi times 2"^ equals the 
number of satisfying assignments of /. Roth [2^ showed that #MON-CNF cannot 
be approximated by a factor of 2'^^ " unless NP = BPP. By above reduction, the 
same inapproximability must hold for the volume of HlLi ^Pi ■ '-' 

This shows that in general there does not exist a polynomial time multiplicative 
e-approximation of the volume of the intersection of bodies in M."^. This holds for 
all classes of objects which include co-boxes, e.g. schlicht domains (cf. Section [273]) . 
Though there is no multiplicative approximation, we can still give an additive ap- 
proximation algorithm, that is, we can efficiently find a number V such that 

Pt[V -eV^ir.<V <V + eV^in] >| 
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Geometric objects 


Volume of the union 


Volume of the 
intersection 


axis-parallel boxes 


#P-liard + FPRAS 


easy 


general boxes 


#P-liard + FPRAS 


#P-liard 


co-boxes 


easy 


#P-hard + APX-hard^ 


schliclit domains 


#P-hard -t- FPRAS^ 


#P-hard + APX-hard^ 


convex bodies 


#P-hard + FPRAS 


#P-hard 



Table 1: Results for the computational complexity of the calculation of the volume of union 
and intersection (asymptotic in the dimension d). 



where V is the exact volume of the intersection and VJnin is the minimal volume of any 
of the given bodies Bi, . . . , Bn- If we could replace Vmin by V in the equation above, 
we would have an FPRAS. This is not possible in general as the ratio of V and Vmin 
can be arbitrarily small. Hence, such an e-approximation is not relative to the exact 
result, but to the volume of some greater body. This is an additive approximation 
since after rcscaling, so that Vnun < 1 we get an additive error of e. Clearly, we get the 
result from above by uniform sampling in the body -Bmin corresponding to the volume 
Knin- Consider V = Knin {Zi + . . . + Zn) /N , where Zi is a random variable valued 1 , if 
the z-th sample point Xi — SAMPLEQuERY(i3min) lies in the intersection of Bi, . . . , _B„, 
and otherwise. Using Chebyshev's inequality one can show quite easily that V gives 
an approximation as desired, if we choose A'' proportional to with the right factor. 
This gives an approximation algorithm with runtime 0(nV{d) + -I- ■^P{d)^, 

yielding C'(^) for (not necessarily axis-parallel) boxes. 



4 Discussion and open problems 

We have proven :^P-hardness for the exact computation of the volume of the union 
of bodies in R'' as long as the class of bodies includes axis-parallel boxes. The same 
holds for the intersection if the class of bodies contains general boxes. We have 
also presented an FPRAS for approximating the volume of the union of bodies that 
allow three very natural oracles. Very recently, there appeared a few deterministic 
pol yno mial-time approximations (FPTAS) for hard counting problems (e.g. 0, 



15l Il5l |29| ) . It seems to be a very interesting open question whether there exists a 
deterministic approximation for the union of some non-trivial class of bodies. Since 
the volume of convex bodies determined by oracles cannot be approximated to within 
a factor that is exponential in d Q , the existence of such a deterministic approximation 
for the union seems implausible. It is also open whether there is a constant C so that 
KMP can be efficiently deterministically approximated within a factor of C, i.e., if 
they are in APX? 

For the intersection we proved that no multiplicative approximation (deterministic 
or randomized) is possible for co-boxes (cf. Definition [T]), but we also presented a very 



^If the integrals are computable in polynomial time (cf. Scction l2.3| l. 

^Theorem |6] even proves that for every fixed e > approximating the volume within 2'^ 
NP-hard. 
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simple additive approximation algorithm for the intersection problem. It would be 
interesting to know if there is a hard class for multiplicative approximation which 
contains only convex bodies. 

Our results are summarized in Table [TJ Note the correspondence between axis- 
parallel boxes and co-boxes. The discrete counterpart to their approximability and in- 
approximability is the approximability of ^DNF and the inapproximability of #SAT. 
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